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ABSTRACT 

Disk fragmentation resulting from the gravitational instability has been proposed as an efficient mechanism 
for forming giant planets. We use the planet Fomalhaut b, the triple-planetary system HR 8799, and the poten- 
tial protoplanet associated with HL Tau to test the viability of this mechanism. We choose the above systems 
since they harbor planets with masses and orbital characteristics favored by the fragmentation mechanism. We 
do not claim that these planets must have formed as the result of fragmentation, rather the reverse: if planets 
can form from disk fragmentation, then these systems are consistent with what we should expect to see. We use 
the orbital characteristics of these recently discovered planets, along with a new technique to more accurately 
determine the disk cooling times, to place both lower and upper limits on the disk surface density — and thus 
mass — required to fo rm these object s by disk fragmentation. Our cooling times are over an order of magnitude 
shorter than those of Rafikov (2005), which makes disk fragmentation more feasible for these objects. We find 
that the required mass interior to the planet's orbital radius is ~ 0.1 M Q for Fomalhaut b, the protoplanet orbit- 
ing HL Tau, and the outermost planet of HR 8799. The two inner planets of HR 8799 probably could not have 
formed in situ by disk fragmentation. 

Subject headings: Instabilities — planetary systems: formation — planetary systems: protoplanetary disks — 
stars: individual (Fomalhaut, HL Tau, HR 8799) 



1. INTRODUCTION 

As more extrasolar planets are discovered, we are increas- 
ingly pressed to describe how planets can form in such a vari- 
ety of environments. Until just recently, observational selec- 
tion biases have resulted in the fact that all observed extrasolar 
planets have been f ound to orbit within a few AU of their star 
dBufler et al.ll2.Q06b . Since it seems unlikely that the se plan- 
ets could have formed in situ (Mayor & Queloz 1995), planet 
migration is usually invoked (lAlibert et all 12005). Unfortu- 
nately, this means that little is known about where — and hence 
how — these planets originally formed. 

In contrast, the technique of direct-imaging has presented 
us w ith a new set of extrasolar planets tha t lie far from their 
star dKalas et alj|2008t iMarois et alj|2008l) . along with a po- 
tential protoplanet dGreaves et al. 2008). Like previous tech- 
niques, direct imaging preferentially detects giant planets of 
several Jupiter masses. Furthermore, planet migration need 
not be invoked to explain how these planets could form at 
their observed locations. 

One possible mechanism for giant plan et formation is core 
accretion follow ed by rapid gas accretion (Pollac k et al J 19961 
llnaba et aT1l2003h . However, this mechanism has difficulty 
forming giant planets at large radii. The primary reason 
for this is that the initial core accretion time scales as r 3 



2009 for recent developments). Provided that the disk sur- 
face density is sufficiently large, this mechanism can form gi- 
ant planetary embryos on time scales of a few orbital periods. 
However, if the surface density is too large, the disk is unable 
to cool sufficie ntly fast for fragmentation to take place at all 
(Rafikov 2005). The combination of these requirements im- 
plies gravitational instability can only form massive planets at 
large radii. 

In this l e tter, we consider the planet Fomalhaut b 
(iKalas et all 1200 8). the triple-planet system HR 8799 
dMarois et al.l 12008}). and the potential protoplanet orbiting 
HL Tau 1 ([Greaves et al. 2008). Each of these systems pos- 
sesses at least one planet with orbital characteristics favored 
by the disk fragmentation mechanism. By determining the 
range of surface densities required to form a giant planet with 
the same semi-major axis as these observed planets, we can 
infer the range of disk masses needed for the fragmentation 
mechanism to have operated in these systems. 

2. DISK FRAGMENTATION 

The s t ability of a thin, massive disk is controlled by the 
iToomrd {T964) Q parameter 



Q = 



where r is the orbital radi us of the planet (llkoma et al]l2000t 
iKenyon & B romley 2008). Thus, while it may take ~ 1 Myr 
to form a gas giant at 5 AU via core accretion, it would take 
~ 1 Gyr for the same process at 50 AU — far longer than 
the o bserved lifetimes of protoplanetary disks (Hais ch et all 
l200l . 

Another mechanism for giant planet formation is disk frag- 
mentation as a consequence of the gra vitational instability 
(lKuipedll95llCameronl|1978blBosslll997L see also the recent 
review by iDurisen et al.ll2007l and IStamatell os & Whitworthl 
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where c s is the isothermal sound speed, is the orbital angu- 
lar frequency (assuming a Keplerian disk), and £ is the sur- 
face density. The disk becomes gravitationally unstable for 
Q < 1. However, even if a disk is gravitationally unstable, it 
can on ly fragment if it possesses a sufficiently short cooling 
time (Ga mmiell2001l:lRice et alJl2003l) . Specifically, fra gmen- 
tation will only occur if 

fcool < Q ! (2) 

1 For our purposes, the distinction between planet and protoplanet is irrel- 
evant, and we will use "planet" in both contexts from here on. 
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where f coo i is the local cooling time for a small, point-source 
perturbation, and £ is a f actor of order unity that can depend 
on the history of the disk dClarke et alj2007l) . We adopt Q < 1 
and £ = 1 for our fragmentation criteria. 

2.1. Local Cooling Tune 

Typically, the effects of cooling have been studied using 
time-dependent hydrodynamic simulations. Inevitably, these 
numerical approaches have to employ significant simplifica- 
tion of the radiation field for the sake of computation time 
(e.g., optically thin cooling or flux limited diffusion). Many 
of the simulations show that frag mentation does occu r given 
sufficiently high surface densiti es dDurisen et alj|2007h . 

In contrast, iRafikovl d2005l) used an analytic, order-of- 
magnitude calculation to show that cooling times derived 
from the equations of radiative transfer were much longer, 
and that fragmentation thus did not work, except at radii 
> 100 AU. Here, we adopt an approach inspired by Rafikov, 
but with a more complete calculation of the radiative transfer. 
In brief, we find cooling times that are, in most cases, over an 
order of magnitude shorter than those given by Rafikov (see 
Nero & Bjorkman, in prep, for a more complete discussion). 
As a consequence, we find that fragmentation over a larger 
range of the outer disk is possible, depending on the details of 
the system. 

We emphasize that the cooling time we calculate here is 
for a perturbation, an d is not th e same as the total disk cool- 
ing time employed by Gammi 1 d200lh . While the later may 
be more convenient for numerical hydrodynamic simulations, 
the former is necessary to properly account for background 
heating by external illumination (i.e., the host star). The per- 
turbation cooling time determines the onset and initial growth 
of the instability (in the linear regime), while the total cool- 
ing time controls the ultimate (typically non-linear) comple- 
tion of the instability. Note, however, that when self-heating 
is small, the perturbation and total cooling times will be the 
same within a factor of order unity. 

The perturbation cooling time f coo i = A£/(87T AHo)> where 
AS is energy per unit area added by the perturbation, and 
A//() is the frequency-integrated Eddington flux at the disk 
surface. We consider an annulus within the disk, which we 
approximate as a plane-parallel atmosphere with finite thick- 
ness. For simplicity, we assume that the perturbation is lo- 
cated at the disk mid-plane and that the disk cools equally 
from its top and bottom surfaces. Under these assumptions, 
the perturbation cooling time is 

-3/4 
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where is the adiabatic constant for the gas, x dlff is the mean 
opacity (absorption plus scattering), r is the optical depth co- 
ordinate, B and AB are the depth-dependent Planck function 
and its perturbation, and B m , T m , and c m are the Planck func- 
tion, the temperature, and the isothermal sound speed at the 
disk mid-plane, respectively. The limits of integration, tq and 
-ro, are the optical depth coordinates at the "top" and "bot- 
tom" surface, respectively. Note that we break convention 
here by placing r = at the disk mid-plane, rather than at 
the top surface. Locations in the disk below the disk mid- 
plane have negative optical depth coordinates, while those 
above have positive. Also, note that we have assumed that 
X dlff is approximately constant over the vertical extent of the 
disk, since the disk is nearly isothermal in the vertical direc- 
tion. While this is not true for the surface layers, the error 



is minimal since most of the disk mass — and thus internal 
energy — is located in the disk interior. Similarly, the verti- 
cally isothermal assumption does not apply when accretion 
is the dominant source of heating (since then there would be 
a significant vertical temperature gradient); however, at large 
radii accretion luminosity is usually not the dominant heating 
mechanism for the disk. 

We assume that the relevant physics (i.e., reprocess- 
ing/absorption of external radiation along with viscous energy 
generation) can be preserved by splitting the intensity into two 
frequency components: diffuse, which corresponds to photons 
that have been reprocessed and emitted by the disk, and exter- 
nal, which corresponds to unabsorbed photons emitted from 
the central star and potentially scattered in the disk. Assuming 
gray opacity (i.e., the appropriate mean for each spectrum), 
the frequency integrated moments of the transfer equations 
are 

dH al K ext rext 

=~^ J (4) 

ext fjext 
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where J is the mean intensity, H is the Eddington flux, and x 
and k are the total and absorptive opacity, respectively. We 
have written all optical depth coordinates in terms of r = 
T m = T ext x diff /x ext - The Eddington factors, / ext = K exl /J ext 
and / dlft = K dlS /J AlS , are assumed to be constant with depth. 
Note that we have included no thermal emission for the exter- 
nal frequency because the disk is typically much cooler than 
the star. 

To find the external radiation, we combine eqs. © and (0, 
which have the solution 



i/ ext = // extii n ML 

sinh/3ro ' 



(8) 



where (3 = K ext ^ ext 



/[(X diff ) 2 / ext ] and H a xt is tne net external 
surface flux. 

There are three sources of energy for the diffuse radiation: 
1) absorption of external radiation -dH ext /dr, 2) accretion 
luminosity L acc with surface flux Hq = (dL- dcc /dA)/8ir, and 
3) the point-source perturbation at the mid-plane AHoS(r). 
Thus, the flux transported by the disk is H = Hq(t/to) + 
A// () sgnT-// ext . 

From eq. (0, we can now obtain the diffuse mean intensity 



J dm = H 
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where g diff = i/ diff (ro)/7 diff (To) is the second Eddington factor. 

Finally, the temperature may now be found from radiative 
equilibrium, eq. (O, which gives 

dh (10) 



B = J dm + X 
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as well as 
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Given the temperature and its perturbation, we calculate the 
cooling time from eq. (O. In terms of the surface density, 
£ = 2to / x dlff , the cooling time 



^cool °C £ 



(12) 



where the exponent a depends on both the disk optical depth 
and the ratio of heating from the star to self-heating from ac- 
cretion (see Nero & Bjorkman, in prep.), but is generally in 
the range 0-2. In the negligible self-heating regime (Hq — » 0, 
B/B m — v 1), eq. © reduces to 
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In the optically thick limit (t > 1 , f m 
ther simplifies to 



1 /3), eq. (Qj) fur- 



^cool ' 



1 



diff ^2 



m lA -\aT*' 



(14) 



Note that this expression is di rectly compara ble to the opti- 
cally thick limit of eq. (4) in iRafikovl d2005l) . with the ex- 
ception that our numerical factor (3/128) is over an order of 
magnitude smaller than his (1/4), thus leading to significantly 
reduced cooling times. This difference from Rafikov's result 
is much larger than the difference between total and perturba- 
tion cooling times. It instead arises from his assumption that 
several factors were of order unity, which is true individually; 
however, the cumulative effect of his assumptions results in 
his making a large overestimate. 

Our analytic expression for the cooling time, eq. ( TBI , as- 
sumes no accretion luminosity and large optical depth. While 
it is probably safe to assume that self-heating is small at large 
disk radii, it is not as certain that the optically thick limit will 
hold. For this reason, we use the full solution of eq. (f3]) when 
calculating t coo i, rather than assuming that this limit holds. 
Nonetheless, eq. (TBl i makes for a good estimate in most cases 
where £ is large enough to satisfy Q < 1 in the outer disk. 

2.2. Fragment Masses 

When determining if disk fragmentation is a viable mecha- 
nism for forming giant planets, another point to consider is the 
issue of producing the proper planetary mass of a few Mj up i te , -. 
While a full treatment of this problem is beyond the scope of 
this letter, we provide a toy model to argue that this is likely 
to be the case. 

If Q < 1 at some radius r, the disk becomes gravitationally 
unstable. Supposing m spiral arms form, each arm has local 
surface density 

w r 

(15) 



£arm = — ~ £ > 

m R 



where £ is the original surface density of the previously ax- 
isymmetric disk, R is the current width of the spiral arms, and 
w is a constant that depends on the winding angle. For sim- 
plicity, we assume that most of the disk mass is confined to 
the spiral arms, while the space between the arms is effec- 
tively empty. In addition, we assume w <~ 1, corresponding to 
moderate winding. 

If fcooi > l/£~2, then the spiral arms are pressure supported 
and are stable against fragmentation. However, if f coo i < 1 /f2, 
then the arms are instead supported by centrifugal forces. 
As R continues to decrease, they will fragment radially once 
the centrifugal support is lost. Balancing self-gravity against 



TABLE 1 

Parameters Used for Each System 



Object 


(Mq) 


R* 
(«©) 


(K) 


r planet(s) 

(AU) 


Fomalhaut b 


2.0 


1.8 


8500 


> 101.5 :l 


HR 8799 b, c, d 


1.5 


1.8 


8200 


68, 38, 24 b 


HL Tau b 


0.3 


0.6 


3700 


65 c 



a IChiang et all <2008T) b IMarois et al] 

1 2008 )° IGreaves et all f2008h 

the centrifugal support, tt 2 R = 7rG£ arm , fragmentation occurs 
when 



R<R f = r 2 J , 



(16) 



where is the stellar mass. The fragment mass is 7r7?y-£ arm , 
so assuming m ~ ir, we find a characteristic fragment mass 



Mf ~ lMjupiter 



lOgcm 



Mr. 



100AU 



(17) 

which is consistent with our requirement to produce Jupiter- 
mass planets. 

3. DISK MASS LIMITS 

The condition that both the Toomre Q and the cooling time 
be sufficiently small (eqs. |Q] and |2)) can be used to place 
limits on what disk surface densities lead to fragmentation. To 
be gravitationally unstable, the Toomre Q condition requires 
that the surface density be larger than a minimum, £ m i n , while 
the cooling condition imposes a maximum surface density, 
£max- Therefore, local disk fragmentation is only possible if 
£min < £ < £max- It °f course follows that for fragmentation 
to be possible at all, £ m j n must be less than £ max . This limits 
the range of radii were fragmentation is even a possibility. 

In Fig.[l] we plot £ m ; n and £ max for the systems Fomalhaut, 
HR 8799, and HL Tau, using the parameters listed in Table[T] 
In all cases we used an accretion rate of 10~ 6 MQyr -1 , a mean 
molecular weight for the disk of p g = 2.33, and an adiabatic 
gas constant of ja = 1 -43. The temperature is determined from 
eq. ( fTOb . using a flared disk model with a power law scale 
height h oc r 5 / 4 , which determines the angle of incidence of 
the ext ernal radiation. We use the dust opacity of Cote ra et"aTI 
(|2001|) and the Rosseland mean gas opacity of Hellin get alJ 
(2000). We note that the stellar parameters in Table[T]are not 
necessarily appropriate if planet formation occurs during the 
Class 0/1 phase when the star is significantly more luminous. 
However, our results are relatively insensitive to this effect, 
since T m is only weakly dependent on the stellar luminosity. 

Each system presented here has a planet that might have 
formed via disk fragmentation, assuming that the local sur- 
face density had a value between £ m ; n and £ max for at least 
a few orbital periods. By assuming a power law for the sur- 
face density, £ oc r~ p , we can calculate a range of disk masses 
(interior to r p i anet ) that satisfies this cond ition. 

A su rvey of 24 circumstellar disks by And rews & Wiliamsl 
(120071) found p « 0.0-1.0 with an average of p ~ 0.5 , while 
the hydrodynamical simulations of Vorobyov & Basul (120091) 
found p w 1 .0-2.0 with an average around p w 1 .5. Disk mass 
limits M c / for p = 0.0, 0.5, 1.0 and 1.5, along with the more 
fundamental surface density limits are given in Table |2] We 
also provide the characteristic fragment mass (approximate 
planet mass) Mf from eq. ( TFTI ) that we would expect from 
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TABLE 2 

Range of Disk Masses that Fragment 



Object 


S 


M d (p = 0.0) 


M d (p = 0.5) 


M rf (p = 1.0) 


M d (p = 1.5) 


Mf 




(g cm -2 ) 


(M S ) 


(M ) 


(M Q ) 


(M ) 


Mjupiter) 


Fomalhaut b 


35-94 


0.13-0.34 


0.17-0.45 


0.26-0.68 


0.51-1.36 


5-21 


HR 8799 b 


62-89 


0.10-0.15 


0.14-0.19 


0.20-0.29 


0.41-0.58 


4-7 


HLTaub 


21-43 


0.03-0.06 


0.04-0.09 


0.06-0.13 


0.13-0.26 


2-4 



< 1 


( 

\ Fomalhaut 




1 


; ' 



1.0 1.5 2.0 2.5 



N 




— 




HR 8799 











1.0 1.5 2.0 2.5 



1 

- N 


1 

HL Tau 


1 





1.0 1.5 2.0 2.5 



logr [AU 



FIG. 1. — Surface density limits for disk fragmentation. The solid line 
denotes E max , which is the maximum surface density for the cooling time 
constraint. The dashed line denotes the minimum surface density for frag- 
mentation, E m j„, which is the locus Toomre Q = 1. Disk fragmentation is 
only allowed in the region S max > £ m j„, which is shaded gray. The locations 
of known planets are plotted as vertical dotted lines. 

the disk fragmentation mechanism. Also note that, in order 
to be conservative, we are using the smallest radius found 
by Chiang et alj (12008b for Fomalhaut b. Using one of their 
better fits (e.g. 115 AU) will decrease our lower disk mass 
limit by a few percent, increase our upper disk mass limit by 
w 30% (which would make fragmentation slightly easier), and 
increase the characteristic fragment mass by w 50%. 

4. DISCUSSION 

We have refined the calculations of Rafikov (2005) and 
found cooling times over an order of magnitude shorter. We 
have used these cooling times, along with the observed stellar 
parameters of Fomalhaut, HR 8799, and HL Tau, to test the 
viability of the disk fragmentation mechanism. We found that 
in each of these systems, at least one planet could have formed 



in situ as the result of fragmentation, assuming the disk mass 
interior to those planets fell within a particular range as indi- 
cated in Table|2] 

While the ranges in Table [2] only span a factor of a few, 
this is not by itself a significant limitation. Even if the local 
surface density is above the upper instability limit, fragmenta- 
tion may still occur since the surface density must eventually 
drop through the unstable regime as the disk evolves and dis- 
sipates. The caveat is that the surface density needs to evolve 
on a timescale longer than an orbital period so that there can 
be sufficient time to fragment. 

Our minimum disk masses for Fomalhaut b, HR 8799 b, 
and HL Tau b are about an ord er of magnitude larger tha n 
those inferred from observations dAndrews & Wiliamsll2007l) . 
Note, however, that this is a problem for all planet forma- 
tion models in general. Even core accretion models require 
an enha nced surface density (although to a s omewhat lesser 
extent) dPollack et alJll996Hlnaba et alj|2003l) . One possible 
mechanism for increasing th e surface density is mass loading 
from an infalling envelope ( Vorobyov & Basu 2006). Con- 
versely, current estimates of disk masses may be too low be- 
cause they depend on: 1) the extrapolation of surface densities 
in the outermost regions of the disk to the inner disk, and 2) 
the rather uncertain dust opacity. For example, larger dust 
grains would require larger disk masses to fit the observed 
SEDs dAndrews & Wiliamsll2007h . 

As further evidence for underestim ated disk masses, numer- 
ical hydrodynamical simulations by Vorobvov (2009) found 
disk m asses much higher than those of I Andrews & W iliams 
d2007l) . In particular, stars like Fomalhaut and HR 8799 can 
support disks as large as 0.5 Mq, while HL Tau could have a 
disk as massive as 0.1 Mq, all of which are within our limits 
for disk fragmentation. We caution, however, that our choice 
of opacity model can have a major effect on our results. For 
example, decreasing the dust opacity raises the temperature 
and decreases the cooling time in the outer disk, resulting in 
disk fragmentation at smaller radii. On the other hand, in- 
creasing the opacity would have the opposite effect. 

Regardless of the above considerations, HR 8799 c and d 
are too close to their parent star to have formed in situ via 
fragmentation under the conditions modeled here. Appeal- 
ing to chronically overestimated dust opacity can only get 
us so far. Dropping the opacity by an order of magnitude 
brings HR 8799 c into the fragmentation zone, but still leaves 
HR 8799 d out. Likewise, twiddling other parameters can 
also move the fragmentation radius inward, but reaching the 
required 24 AU with reasonable parameters does not seem 
possible. We therefore conclude that HR 8799 c and d likely 
did not form in situ as the result of disk fragmentation (of 
course, they could have formed at larger radii were the disk is 
more likely to fragment and migrated inward). 

For those planets that could form by disk fragmentation, 
we find characteristic fragment masses (approximate planet 
masses) of a few Mj up i tel for the lower end of unstable disk 
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surface density. Our estimates are mostl y consistent with ex- 
pectations, although Chi ang et al.l (120081) found Fomalhaut b 
to have an upper mass limit of 3 Mj up i tel -, which is 60% lower 
than our lowest estimate of 5 Mj up i tel . Nonetheless, consider- 
ing the crudeness of eq. ([T7l >. we are not convinced that this 
discrepancy rules out Fomalhaut b from having formed as a 
result of disk fragmentation. 

In conclusion, given our current uncertainty of typical disk 
masses and dust opacities, we cannot rule out planet forma- 
tion from disk fragmentation. Of the three currently known 
systems that might posses planets formed in this manner, all 



are able to produce at least their outermost giant planet via 
fragmentation, given a large enough disk mass at some point 
in their history. 
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